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Abstract 

In this paper, we determine the density functions of doubly noncentral singular 
matrix variate beta type I and II distributions. 

1 Introduction 

Matrix variate beta type I and II distributions h ave been studied b y many authors us - 
ing different definitions i n the nonsingular case, se e Kshirsagail ( 1961 ). Constantine ( 1963 ). 

Olkin and Ru binMl964l). |jamed (ll963).ISrivastaval (ll968h.lKhatril (ll970l ). ISr^ and Khatri 

( 1979 ). Cadet (1990} and Diaz-Garcia and Gutierrez- Jaimed ( 2007 ) among ma n y others. 

Recen tly, t hese distributions have been studied in the s ingular cases by lUhligj (|l994l ). lDiaz-Garda and Gutierrez 



( 19971 ) and lDiaz-Garcia and Gutierrez- Jaimed ( 2008a ). Beta type I and II distributions play 
a very important role in several areas of multivariate statistics, such as ca nonical correlatio n 
analysis, the gen eral li near hypothesis in MANQV A and shape theory, see iMuir head (Il982l) . 
Srivastaval (|l968[) and lGoodall and Mardial (|l993l) . 

In the nonsingular case in particular in the nonsingular case, the doubly noncentral 
density functions of matrix variate beta type I and II distributions have been studied by 
diverse authors, but with special emphasis on symmetrised density functions and their 
application to the theory of matrix variate distribut ion, the multivariate Behrens-Fishe r 
problem and the generalis ed re gression coefficient, see Chikuse (1980. 198ll). David ( 1980l ). 
Chikuse and Davis ( 19861 ) and Diaz-Garcia and Gutierrez-Jaimea (■2008a): ?. 

Using Greenacre's definition of the symmetrised density function ( Greenacr3 . 1973f) . in 
an inverse way, we obtain the doubly noncentral nonsymmetrised density functions or 
simply, the doubly noncentral density functions of the singular matrix variate beta type I 
and II distributions, see Section [3l Moreover, as particular cases we find the noncentral 
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density function of the singular matrix variate b eta type I and II dis tributions, fro m where 
we re solve, indirectly, the integral prop osed by IConstantind (Il963l). discussed by Khatri 

1970 ) and reconsidered in (|Farrel]| . ll985l p. 191), see also lDiaz-Garcia and Gutierrez- Jaimez 

20071 ). in singular and nonsingular cases. 



2 Preliminary results 

In this section we give some definitions and notations for the singular matrix variate beta 
type I and II distribution. We also include two results for the symmetrised function and 
invariant polynomials with matrix arguments. 



2.1 Singular beta distributions 

Consider the following definition and notation. 

Let C be a non-negative definite mxm; then C^/^(C^/^)' = C is a reasonable nonsingular 
factorization of C, and in particular C^^^ can be m x m upper-t r iangular matrix or an 
mxm non- negat ive definite square root, see Gupta and Nagar ( 2000l ). Srivastava and Khatril 
(|l979t ) and .MuirheadI (Il982l) . 

Let A be an m X m non-negative definite random matrix with Pseudo-Wishart distri- 
bution with r degrees of freedom and a symmetric matrix of parameters S. We then state 
that A ~ VWm{r, S), Re(r) < (to - 1). If Rejr) > (to - 1) then A is said to ha ve a 
Wishart distribution, with A~ >Vm(s, S), see (|Muirheadl . Il98l p. 82), lUhligl (|l994l ) and 
Diaz-Garcfa and Gutierrej ( 1997? ). 



Definition 2.1. Let A and B be independent, where A ^ 'pyVm(r, I) and B ~ yVm{s,T). 
We def ine U = (A-^B)-l/^A((A + B)-l/ ^)^ Then its density function is given and denoted 
as (see Diaz-Garcia and Gutierrej ( 19971) 1 



S/,„(U;(7,r/2,s/2) = c|L 



(r-m-l)/2 



(s-m-l)/2 



(dU), 0<U<I„ 



(1) 



U is said to have a singular matrix variate beta type I distribution, and this is denoted as 
U - S/,„(g,r/2, s/2), Re(s) > (to- 1); where U = HiLH;, with Hi e V,,™; V,,™ = {Hi S 
sfjmxg|jj/^jj^ = IJ denotes the Stiefel manifold; L = diag(;i, . . . , Ig), I > h > ■ ■ ■ > Iq > 0; 
q — m (nonsingular case) or g = r < to (singular case); 



7r(-""-+'-g)/2r^[(r + s)/2] 
r,[r/2]r™[s/2] 



(2) 



(dU) denotes the Hausdorff measure on {mq — q{q — l)/2)-dimensional manifold of rank-q 
positi v e sem idefi nite mxm matrices U with q dis tinct nonnuU eigenvalues, given by (see 
lUhliel (|1994| ) and iDiaz-Garcia and Gutierrez! (|l997l )) 



(dU) = 2-'n/r"'n(^^ - A * ^ (HWHi), 



(3) 



i<j 



where (H'ldHi) denotes the invariant measure on Vq^m and where r,„[a] denotes the multi- 
variate gamma function, this being defined as 

r™[a] = / etr(-R)|Rr-("+i)/2(dR), 

"'R>0 

Re(a) > (to - l)/2 and etr(-) = exp(tr(-)). 
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Similarly, we have: 

Definition 2.2. Let A ^ PWm(r, I) and B ~ Wm(s,I) be independent. Furthermore, 
let F — B^^/^A(B^^/^)'. F is then said to have a singular matrix variate beta type II 
distribution, denoted by F ~ BII„i{q,r /2, s /2) and if F = HiGH'^, with Hi G V,,™ and 
G = diag(ffi, . . . ,,(7o); gi > ■ ■ ■ > gq > 0, its density function is given and denoted as (see 
Diaz-Garcia and Gutierred (|l997l) ) 



i3//„(F;q,r/2,s/2) = c|G|('-""i)/2|j^^P|-(r+.)/2(^p)^ p > 0^ 

where c is given by ([2]), Re(s) > (m — 1) and (dF) is given in an analogous form to ([3]). 

Let us now extend these ideas to the doubly noncentral case, i.e. when A ^ 7^W„i(r, I, Oi) 
and B ~ Wm(s, I, J^2)- In other words, A has a noncentral Pseudo-Wishart distribution 
with a matrix of noncentrality parameters rii and B has a noncentral Wishart distribution 



with a matrix of noncentrality par a meters 112, see lDiaz-Garcia et al\ (jlQQTT) . Subsequently, 



Diaz-Garcia and Gutierrez- Jaimez ( 2008al ) reported the following 



Lemma 2.1. Suppose that U has a doubly noncentral matrix singular variate beta type I, 
which is d enoted as U ^ BI„i{q, r/2, s/2, Oi, 172). Then using the notation for the operator 
sum as in lDavi^ ][198(\ } its symmetrised density function is found to be 

dF,{lJ) = BImiV; q, r/2, s/2) etr (-^(Oi + n^)) 

^ + ^)), ^»2)c;-^(u, (I - u)) 

TOi/i < U < I, Rc(s) > (m — 1), (a)r is t/ie generalised hypergeometric coefficient or 
product of Pochhammer symbols and Cl'^j- , ■ ) denotes the i nvar iant polynomials with matr ix 
arguments defined in Mavis (1980. ), see also lChikuse i '198(\ ) and Chikuse and Davis 



Moreover: 

Lemma 2.2. Suppose that F > has a doubly noncentral singular matrix variate beta 
type II, which is denoted as F r/2, s/2, r2i, r22). Then its symmetrised density 

function is 

dG,{¥) = BII„,{¥- q, r/2, s/2) etr {-^{n, + n^)) 

^ ^(r + s), ^02)C;-^((I + F)-iF, (I + F)-i) ^^^^ 

where F > and Re(s) > {m — 1). 

2.2 Symmetrised function and invariant polynomials with matrix 
arguments 

Consi der the follow extension of the definition given by iGreenacrd ( 19731 ). see also Roux 



(119751 ): 



Definition 2.3. The symmetrised density function of the non-negative definite matrix X : 
m X m, which has a density function f^{lC), is defined as 

/,(X)=/ /,(HXH')(dH), HeO(m) (5) 

Jo{m) 
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where 0{m) = {H G 5R'"^'"|(HH') = (H'H) = I,„} and (dH) denotes the normaHsed 
invariant measure on 0(to), then 

/ m) = 1 

JO(m) 



(|Muirheadl . Il982l . p. 72). 

Now consider the foll owing theorem , whic h generalises eq. (5.4) of David ( 1980( ). and 
proof of which is given bv lDfaz-Garcfa I (|2006l) . 



Lemma 2.3. Let A, B, X and Y be m x m symmetric matrices, then we have 



I C^'^(AH'XH,BH'YH)(dH) 

Jo(m) 



c;'^(A,B)c;'-^(x,Y) 



K,\j 



(6) 



with 



3 Doubly noncentral singular matrix variate beta dis- 
tributions 

Taking into account equation ^ it is now is possible to propose an expression for the 
(nonsymmetrised) density functions of doubly noncent ral matri x varia t e bet a type I and 
II distributions, applying the idea of iGreenacrd (|l973l ) (see alsolRoujJ (119751 )1. but in an 
inverse way. 

Theorem 3.1. Assume that U ^ BIm{q, r /2, s/2, fJi, ^2)- Then its density Junction is 
dF^ (U) = BlmiV; q, r/2, s/2) etr (-i(f2i + f^a)) 

X E ft, n7,;„ gr(l"iU,ln2(I~U))(dU) 



with < U < I, Re(s) > (m - 1). 



Proof. First observe that BIm{V; q, r/2, s/2) is a symmetric function, then ;B/m(HUH'; q, r/2, s/2) 
BI,^[\J-q,r/2,s/2). Thus 



dF,{\]) = 6/,„(U; q, r/2, s/2) etr (-i(Oi + il^)) 



for a function h. By ([6]) observe that 



/i(HUH')(dH)(dU), (7) 



0(m,) 



( 1 



^l.HUH', - HUH')) (.H) ^ ^^(M^ d " U)) 



^C,(I) 



Then, by applying ([5]) in an inverse way, in ([7]) we have 

/i(U) = 0^'^ C^'^ (iJ^iU, if22(I - U)) 
from where the desired result is obtained. 



□ 
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Theorem 3.2. Assume that F ~ BIIm{q,r /2,s/2,D,i,D,2). Then we find that its density 
function is 

dG^(F) = BII^{¥- q, r/2, s/2) etr + ^2)) 

X E yfW^^r G^i(i + F)-F,ir2.(i + F)-)(dF), 

where F > aKC? Re(s) > {m — 1). 

Proof. The proof is paraUcl to that given for Theorem 13.11 □ 

In the foUowing two coroUaries we shaU obtain, as particular cases, the noncentral density 
functions of singular matrix variate beta type 1(A), 1(B), 11(A) and 11(B) distributions. 

Corollary 3.1. Under the hypothesis of Theorem lS.ll 

i) // rii = 0, i.e. A ^ 'PWm{r,T), then we obtain the noncentral singular matrix variate 
beta type 1(A) distribution denoted as 

V ^BI{AU{q,r /2, s/2, 
Its density function is then given by 
dFuiU) = i3/„,(U; q, r/2, s/2) etr {-^n^) 

X iFi (i(r + s); is; iJl2(I " U)) (dU) 



ii) Alternatively, if ^2 — 0, i.e. B ~ VVm(s,I), then we obtain the noncentral singular 
matrix variate beta type 1(B) distribution denoted as U ^ BI{B)m{q, r/2, s/2, ^i), for 
which its density function is 

dF^{V) = Bl^iV; q, r/2, s/2) etr {~^n,) (i(r + s); is; i^iU) (dU) 

with < U < I, Re(s ) > (m — 1) an d where i^i(-) is the hypergeometric function with 
matrix arguments, see iMuirhea^ . 198^ . definitions 7.3.1, p. 258). 

Proof. The density functio ns in two item s are a consequence of the basic pro perties of 
invar iant polynomials, see ( Davis . 19791 equations (2.1) and (2.3)), see also ( Chikuse 



1980l equations (3.3) and (3.6)). □ 



Corollary 3.2. Under the conditions of Theorem \3.2[ 

i) if ft I = 0, i.e. A ~ VWmir,!), then we obtain the noncentral singular matrix variate 
beta type 11(A) distribution denoted as 

F ^ BI I iA)raiq, r/2, 3/2,^2), 

and its density function is 

dGF{F) = BIIraiF; q, r/2, s/2) etr (-^^2) 

X iFi (i(r + s); is; 1^2(1 + F)-i) (dF) 
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ii) iffl2 — 0, i-e. B ~ VVm(s,I), then we obtain the noncentral singular matrix variate beta 
type 11 (B) distribution denoted as 

F^S//i(S)„(g,r/2,s/2,f2i), 

for which its density function is 

dG^ (F) - BIIraiF; q, r/2, s/2) etr (-iJ^i) 

X iFi (i(r + s); is; + F)-1F) (dF) 

with < F and Re(s > (to - 1)). 

Proof. The proof is analogous to that given for Corollary 13. II □ 

Conclusions 



IChikusel (|l98Ci ). [Chikuse (1981) and Davis ( 19791 ) have found the symmetrised doubly non- 



central density functions of the nonsingular matrix variate beta type I and II distributions. 
However, the question of nonsymmetrised density functions (or simply density functions) 
remained to be resolved. In t his paper, by applying Greenacre's definition of symmetrised 



function of I Greenacrd . Il973l) in an inverse way, we respond to these two open problems 



with res pect to singular a nd nonsingular cases. Furthermore, in another w ay to the method 
given inlPiaz-Ga rci a and Gutierrez- Jaimez (2007) and in Dia z-Garcia and G utierrez- Jaimcj 
(|2008ah . we obtain the noncentral density functions of singular matrix variate beta type 



1(A), 1(B) , 11(A) and 11(B) distributions, from where, implicitly, we resolve the integral pro- 
pos ed bv IConstant inc (1 963), Khatr i (1970) and reconsidered in (Farrell, 1985, p. 191), see 
also iDiaz-Garcia and Gutierrez- Jaime j (|2007l ). in the singular case, and in the nonsingular 
one, of course, by simply taking q = m. 
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